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ABSTRACT 

Aims. NGC 2264 is a young cluster with a rich circumstellar disk population which makes it an ideal target for studying the evolution 
of stellar clusters. Our goal is to study the star formation history of NGC 2264 and to analyse the primordial disk evolution of its 
members. 

Methods. The study presented is based on data obtained with the Infrared Array Camera (IRAC) and the Multiband Imaging 
Photometer for Spitzer (MIPS) on board the Spitzer Space Telescope, combined with deep near-infrared (NIR) ground-based 
FLAMINGOS imaging and previously published optical data. 

Results. We build NIR dust extinction maps of the molecular cloud associated with the cluster, and deteiinine it to have a mass of 
2.1x10^ Mq above an Ay of 7 mag. Using a differential Kj-band luminosity function (KLF) of the cluster, we estimate the size of 
the population of NGC 2264, within the area observed by FLAMINGOS, to be 1436±242 members. The star formation efficiency is 
>~25%. We identify the disk population and divide it into 3 groups based on their spectral energy distribution slopes from 3.6^m 
to ipm and on the 24 /jm excess emission: (i) optically thick inner disks, (ii) anaemic inner disks, and (iii) disks with inner holes, 
or transition disks. We analyse the spatial distribution of these sources and find that sources with thick disks segregate into sub- 
clusterings, whereas sources with anaemic disks do not. Furthermore, sources with anaemic disks are found to be unembedded (i.e., 
with Ay <3 mag), whereas the clustered sources with thick disks are still embedded within the parental cloud. 

Conclusions. NGC 2264 has undergone more than one star-forming event, where the anaemic and extincted thick disk population ap- 
pear to have formed in separate episodes: the sources with anaemic disks are more evolved and have had time to disperse and populate 
a halo of the cluster We also find tentative evidence of triggered star-formation in the Fox Fur Nebula. In terms of disk evolution, 
our findings support the emerging disk evolution paradigm of two distinct evolutionary paths for primordial optically thick disks: a 
homologous one where the disk emission decreases uniformly at NIR and mid-infrared (MIR) wavelengths, and a radially differential 
one where the emission from the iimer region of the disk decreases more rapidly than from the outer region (forming transition disks). 

Key words, stars: circumstellar matter - stars: formation - ism: clouds 



1. Introduction 



The Spitzer Space Telescope JWerner et al.l l2004h . a space- 
borne observatory of unprecedented sensitivity and spa- 
tial resolution in the infrared, has brought forth an excit- 
ing new era in the study of star formation. In particu- 
lar, the Multiban d Imaging Photometer for Spitzer (MIPS) 
iRieke et al.' '2004") has excelled in identifying protostars (e.g. 
Muzerolle et al. 2004; Teixeira et al. 2006; Muench et al. 2007; 



Winston et aP l2007h . and the Infrared Array Camera (IRAC) 
( Fazio et al.ll2004b|) has permitted the detailed s tudy of jets (e.g. 



iNoriega-Crespo et al.l2004t Teixeira et al.l2008l) and circumstel- 



lar disks (e.g. Sicilia-Aauilar et al. 2006; Hernandez et al. 2008; 
[Luhman et al. 2008 ; Sung et al. 2009 ; Sung & Bessell 2010 ). 
One of the earliest Spitzer studies of circumstellar disks was 
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that of iLada et al.l (l2006h . who conducted a detailed study of 
the frequency and nature of circumstellar disks in IC348. The 
study found that this 2-3 Myr old cluster has a disk fraction 
of 50%+6%, and that this frequency is a function of spectral 
type, where K6-M2 stars possess the greatest disk longevity. 
Regarding the nature of the disks, iLada et all (l2006h identified 
optically thick and optically thin or anaemic inner disks accord- 
ing to their spectral energy distribution (SED) slope from 3.6//m 
to 8/im (aiRAc)- The majority of the IC 348 sources classified as 
classical T Tauri stars (CTTS) (i.e., accreting according to the 
equivalent width of their Ha emission) were found to possess 
thick disks, while the majority of the weak-line T Tauri stars 
(WTTS) (i.e., non-accreting according to the equivalent width 
of their Ha emission) were found to be disk-less, leading the au- 
thors to conclude that the gaseous and dust components of disks 

evolve on similar timescales. 

NGC 2264 is a young hierarch ical cluster (iPichel Il993l 
lLadaetal.iri993l; iLada & Lad"all2003h located in the Mo noceros 
OBI association, that has been extensively studied since I Walked 
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(Il956l) 's seminal work identifiying its pre-main sequence pop- 
ulation (PMS). The most recent distance measurement to the 
cluster places it at 913 +110pc (Baxter et al. 2009) (the value 
used in this paper); previous estimates of the cluster distance had 
varied between 700 pc a nd 800 pc. The oldest members of the 
cluster are ~5Myr old (iDahm & Simonll2005h , although there 
is significant observational evidence of ongoing star formation 
such as CO outflows (e.g. Margulis et al. 19881), Herbig-Haro 
objects ( iReipurth et al.ll2004l) . and luminous far-i nfrared sources 
detec ted by the Infrared Astronomical Satellite (iMargulis et alJ 
Il989l) . The PMS po pulation has furtherrn ore been associated 
with X-ray emission (Flaccomio et al. 2006) an d active accretio n 
(Ha observations, Dahm & Simon 2005: Furesz et al.l l2006h . 
NGC 2264 is one of the classical target clusters that has time 
and time again given valuable insights into the star formation 
process of low- and intermediate-mass stars. Initial results ob- 
tained from Spitzer observatio ns of NGC 22 64 led to the dis- 
covery of the Spokes cluster (iTeixeira et al.l 12006.) , where the 
nearest neighbour separations between its class I members are 
consistent with thermal Jeans fragmentation of the parental fila- 
mentary cloud. Within the Spok es cluster, a dense micro-cluster 
of protostars was also revealed (lYoung et al.ll2006]). and follow- 
up ob servations with the Submillimeter Array (ITeixeira et alj 
|2007|) found new Class sources whose average separations 
were also consistent with the Jeans length for the parent molec- 
ular core^ffli£micr{>cluster More recently, jSung et al. (2009) 
and lSung & Besselll ( |20IO|) reported the spatial distributions of 
Class I vs. Class II sources; they confirmed the spatial con- 
centration of Class I sources in the Spokes cluster and identi- 
fied another grouping of Class I sources near IRS 1, that coin- 
cides with a c luster i ng of compact millimetre sources found by 
iPeretto et al.] (l2006i l2007h . The northern region of NGC 2264, 
near the 07 star S Mon, is mostly compr i sed of Class II sources; 
ISung et al.1 (l2009l) and ISung & Bessdll (l2010l) found that the 
fraction of primordial disks increases with distance from S Mon, 
indicating that UV radiation plays an important role in disk evo- 
lution. The cluster has also been observed by the space-borne 
telesco pe COROT, which p rovided exquisitely detailed light- 
curves dAlencar et al.l 120101) . Comparison of the three different 
types of identified light curves (AA Tau-like, spot-like, and ir- 
regular) with Spitzer IRAC d ata showed remarkable correlation 
and led lAlencar et aP (l2010l) to conclude that -30 to 40% of 
sources with inner disks are associated with AA Tau-like light 
curves. In this paper, we revisit the NGC 2264 cluster by char- 
acterising the disk population and correlating its properties with 
those of the molecular cloud from which it formed. 



2. Data 

2.1. Spitzer IRAC and MIPS data 

For a description of the acquisition of the IRAC data u s ed, w e 
refer the reader to ITeixeira et al.l(l2006l) and lYoung et al.l(l2006h. 
The source extractor SExtractor software (iBertin & ArnoutsI 
Il996l) was used to generate a list of point sources for each wave- 
length. Aperture photometry was performed on these detected 
sources using the Image Reduction and Analysis Facility (IRAF) 
routine APPHOT, with an aperture radius of 2 re-sampled IRAC 
pixels, and a sky annulus of inner and outer radii of 2 and 6 
re-sampled IRAC pixels, respectively. The zero point magni- 
tudes used for the IRAC pass-bands of effective wavelengths 
3.6fim, 4.5 ;um, 5.8/Ym, and 8/Ym were 17.83, 17.34, 16.75, and 
15.93 mag, respectively. IRAC band merging was performed us- 



ing the IDL function srcor Q, where we used a matching radius 
of 1 pixel (0.9")- The photometry of bright saturated sources 
extracted from the 10.4 s integration images were replaced by 
the photometry obtained from the short-exposure, 0.4 s integra- 
tion images. The IRAC data used in this paper was limited to 
photometry with errors smaller than 0.1 mag. The completeness 
limits are 13.25 mag for the 3.6 fim, 4.5/^m, and 5.8yum IRAC 
passbands, and 12.75 mag for 8/Ym IRAC passband. The errors 
are systematically larger for the 5.8 //m and, especially, the 8 fixn 
band photometry. For a description of the MIPS data used (ac- 
quisition, calibration , an d photometry), we refer the reader to 
ITeixeira et al.l (2006) and lYoung et al.l(l2006l) . The completeness 
limit of the MIPS 2A jjm. data is 7.0 mag; we do not use the 
70yum or 160jum data in this paper as very few pre-main se- 
quence sources are detected at these longer wavelengths. We 
merged these data with the IRAC sample using the same method 
described above for the IRAC band-merging, using the closest 
match within a radius of 3.5". 



2.2. FLAMINGOS imaging 

Deep near-infrared (NIR) / (1.25yum), H (1.65//m), and 
K (2.2 yum) data was obtained with the FLoridA Multi- 
object Imaging Near-IR Grism Observational Spectrometer 
(FLAMINGO^ lEIstonll 19981: iLevine et al.ll2006h mounted on 
the Kitt Peak National Observatory (NOAO) 2.1m. The data 
is part of the FLAMINGOS giant molecular cloud survejy. 
FLAMINGOS employs a 2048x2048 HgCdTe HAWAII-2 in- 
frared array with 18 ;um pixels. This corresponds to a plate scale 
of 0.608" pixel' and a field of view of 20.5'x20.5'. Table [U 
summarises the observing information for the fields acquired 
in NGC 2264. The FLAMINGOS observations cover an area 
of w 0.3 square degrees. The seeing varies within the dataset, 
at all wavelengths, from 0.9 to a;1.5" FWHM. The data re- 
duction, calibration, and photometry was performed with cus- 
tom FLAMINGOS data analysis and reduction pipelines within 
IRAF. For a detailed descrip tion of these, we refer the reader to 
iRoman-Zuniga et al.1 (|2008|) . 

We compared the FLAMINGOS and Two Micron All Sky 
Survey (2MASS) (iKleinmann et alJll994l:ISkrutskie et al.ll2006i) 
data and identified the sources with saturated FLAMINGOS 
photometry (the bright sources have FLAMINGOS photometry 
that is systematically fainter then their corresponding 2MASS 
photometry). These saturated sources had their FLAMINGOS 
photometry replaced by 2MASS photometry. The saturation lim- 
its for FLAMINGOS derived from this comparison are 11.0, 
1 1 .5, and 12.0 mag for J, H, and K, respectively. 

The area covered in NGC 2264 by the FLAMINGOS survey 
is smaller than that covered by the IRAC observations. We there- 
fore made use of 2MASS photometry to complete the NIR cov- 
erage of the region observed by IRAC. The NIR data are thus not 
uniform in depth since the FLAMINGOS data has completeness 
limits of 19.0, 18.5, and 17.5 mag for /, H, and K,, respectively, 
while the shallower 2MASS data have corresponding complete- 
ness limits of 15.5, 14.5, and 13.5 mag. The analysis presented 
in this paper is based on NIR photometry that is within the com- 
pleteness limits and has photometric errors smaller than 0. 1 mag. 
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from software from the Ultraviolet Imaging Telescope. 

^ see http : //flamingos . astro . uf 1 . edu 

^ NOAO Survey Program 2000B-0028: Towards a Complete Near- 
Infrared Spectroscopic Survey of Giant Molecular Clouds (PI: E. Lada) 
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2.3. Ancillary data 

To perform a more complete study of NGC 2264, we used pre- 
viously published datasets, in particular the ancillary data cat- 
alogue (ADC) built by Rebull et al. (2006), which is an ex- 
tensive collection of multiwavelength photometric and spectro- 
scopic data of NGC 2264. 

We added additional photometric and spectroscopic datasets 
to the ADC. One of the additional datasets comprises the data 
on NGC 2264 collected bv Mevers-Rice ( 1995) ; these consist of 
VRI photometry and optical spectroscopy. Spectra of 361 stars 
were obtained in January 1993 using the Boiler & Chivens spec- 
trograph (3700 A to 7000 A) with the Steward Observatory 2.3 m 
Bok telescope at Kitt Peak, Arizona. A first set of imaging data 
were acquired during November 1989 and January 1990 using 
the Bok telescope and a 800 x 800 pixel CCD, and a second set 
was obtained in December 1992 with the same telescope but us- 
ing the Loral/Fairchild camera (2048 x2048; the data obtained 
in the second set was binned to match the plate scale of the 
first set). The photometric data were calibrated using data for 
the globular clusters NGC 7790 and NGC 2419. 

An X-ray study of NGC 2264 was conducted by 
iFlaccomio et al.l (l2006l) using ACIS-I onboard the Chandra 
X-ray Observatory. The data obtained has a total integration 
time of 97 ks (observation ID 2540; general observer proposal 
PI S. Sciortino), and the area observed covers 17' x 17' centred 
on (Q',5)(y2000)=(06''40"58.7\-H09°34'14") resulting in the 
detection of 420 sources. These data were subsequently added 
to the ADC. 

The fina l datas et added corresponds to that published by 
iFtiresz et al.l (|2006|) . These authors performed a radial veloc- 
ity survey of NGC 2264 by obtaining spectroscopic observa- 
tions of 990 optically visible stars selected from the studies de- 
scribed above. The data were obtained using the Hectochelle, 
mounted on the 6.5 m MMT in Mount Hopkins, Arizona 
(U.S.A). [Furesz et al. (2006) classified 471 stars as cluster mem- 
bers (according to criteria based on their radial velocities and/or 
Ha line profiles), of which 344 have well-determined radial ve- 
locities. 

The analysis we present in this paper was done using a subset 
of this final merged catalogue: only sources that were detected in 
all four IRAC bands were considered, and we refer to this subset 
henceforth as the IRAC sample. 



3. Results 

3.1. Identification of the disk population 

To begin characterising the disk populatio n of NGC 2264, we 
uses a tool developed by iLada et aD (l2006l) in their study of the 
young cluster IC348; they measured the slope of the SED, for 
each source, between 3.6/im and 8;um. We used this same mea- 
sure, aiRAC, to classify the evolutionary status of sources with 
respect to their inner circumstellar material (emission longwards 
of 8 /zm is considered to arise from outer circumstellar material) . 
The sources in the IRAC sample were classified according to 
the aiRAC empirical scheme shown in Table [2] Cl ass I sources 
viewed pole-on may have flat spectrum SEDs dCalvet et al.l 
119941: JRobitaille et alJl2006l) . and so may sources with optically 
thick circ umstellar disks (CTT S or Class II sources) if viewed 
edge-on (JRobitaille et al.l 120061) . We thus added an intermedi- 
ate evolutionary phase between Class I sources and s ources with 
thick d isks to the empirical classification scheme of iLada et al.l 
(l2006l) : the colloquially termed ^ar spectrum phase. Several au- 



thors have suggested that flat spectrum sources may indeed be a 
distinct evolutionary group altogether, such as JGreene & Ladal 
(l2002b who found that the accretion rates of flat spectrum 
sources are consistently higher than those for CTTS and lower 
than those for Class I sources, and iMuench et al.l (|2007|) . who 
found that flat spectrum sources are intrinsically more luminous 
than Class II sources. 

Figure[T]shows the distribution of q-irac for 1437 IRAC sam- 
ple sources (open hist ogram) and for 349 candidate members 
previously identified bv lFuresz et al.l (l2006l) from their radial ve- 
locity survey (filled hatched-line histogram). The IRAC sample 
ffiRAC histogram shows that the majority (985) of the sources 
have naked photospheres - these are comprised of disk-less PMS 
sources of NGC 2264, as well as foreground and background 
field stars. The IRAC sample has 408 sources that have circum- 
stellar material (Q'T RAr> -2.56), 160 of which had been previ- 
ously identified by iFtiresz et al.l (|2006|) as candidate members. 
Of the 349 sources, »46% have circumstellar material, either 
fr om envelopes (there are 1 1 Class I and flat spectrum sources 
in IFtiresz et aLJ (l2006l) 's list that are detected in all four IRAC 
bands), thick (30%), or anaemic disks (12%). We note that the 
statistics for flat spectrum and Class I sources are particularly 
incomplete since these are very red sources and some may not 
have been detected in all four IRAC bands; the above fraction of 
sources with circumstellar material is thus an underestimate. 

Extragalactic sources have infrared excess emission that can 
be erroneously inter preted as a signature o f a star with a cir- 
cumstellar disk (e.g. |j0rgensen et al.l l2006h . Although the ex- 
tensi ve molecular clou d associated with the NGC 2264 cluster 
(e.g. 'Ri dge et al.ll2003l) is relatively efficient in screening against 
the contamination of extragalactic sources in the IRAC sam- 
ple, it is still necessary to prune the dataset to reduce the num- 
ber of these interlopers. iFazio et aP (l2004al) found that the con- 
tamination by extragalactic sources, such as PAH-rich galax- 
ies, is greater than 50% for sources fainter than 14.5 mag at 
3.6/zm. To reduce the contamination of extragalactic sources, 
we removed 8 sources from our IRAC sample that are fainter 
than 14.5 mag at 3.6/zm. This extragalactic decontamination 
technique has also been used b y other authors to filter their 
data (e.g. iHemandez et al.l 120071 in their study of the young cr 
Orionis cluster). Figure |2] shows that the luminosity functions 
for V- and J-bands are similar for sources with thick disks (open 
histogram), and sources with anaemic disks (filled histogram). 
Since the thick and anaemic disk populations have consistent 
luminosity functions, we conclude that the contamination of ex- 
tragalactic sources in the anaemic disk population is very low. 
Futhermore, we found that the average radial velocity of sources 
with anaemic disks is 18.9+2.8kms"', while sources with thick 
disks have an average radial velocity of 20.0+4.6 km s"'. The 
consistency of the average radial velocities of the disk popula- 
tions also led us to conclude that there is little contamination by 
foreground and background stars. 

The cleaned IRAC sample, used henceforth, thus has 1401 
sources of which 12 are Class I sources, 28 are flat spectrum 
sources, 214 are sources with thick disks, and 154 are classified 
as sources with anaemic disks. Disk-less sources that are heavily 
reddened have ffiRAC values that would lead to t heir identifica- 
tion as sources with anaemic disks, as shown by iMuench et alJ 
(l2007l) . The authors used the photosphere of a KO diskless star 
and created reddened SEDs for various amounts of extinction; 
they found that for Ay > 10 mag, the SED slope q-irac of a back- 
ground diskless star was greater than -2.56, and thus the red- 
dened star could be erroneously classified as having an anaemic 
inner disk. Indeed 36 such sources were found in NGC 2264, so 
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the number of sources with anaemic disks is actually considered 
to be 118. 

The IRAC sample has 386 sources with previously deter- 
mined spectral types. Figure |3] shows how aiRAc varies as a 
function of spectral type. As the diagram shows, the majority 
of sources with thick or anaemic disks are late-type (>G9), al- 
though some early-type sources have thick disks. A comparison 
of Figure [3] with Figure 1 of Lada et al. ( 2006) shows great sim- 
ilarity in the variation of ffiRAC with spectral type between the 
NGC 2264 and IC 348 clusters. If these sources are coeval, this 
diagram indicates that the dis ks around lower m ass stars may be 
longer lived, as suggested bv lLada et al.l (l2006l) . 

3.2. The molecular cloud 

To better characterize the nature of the disk population of 
NGC 2264, it is useful to study the structure of the molecular 
cloud associated with the young cluster We began this study by 
building a NIR dust extinction map of the region surrounding the 
NGC2264 cluster, shown in Figure |4] The map was built using 
the near-infrared colour excess (NICE) method (see Lada et al.l 
II994tlAlves et al J 19981 for a detailed description of the method) 
and 2MASS H and K, data dSkrutskie et al.l2006h . Basically, the 
NICE method measures the reddening of the (H-Ks) colour of a 
star background to a molecular cloud. Young stars with circum- 
stellar disks may also have reddened (H-Ks) colours due to the 
reprocessed stellar emission from the dust in their inner disks - 
it is therefore necessary to remove these disked sources from the 
NIR photometric sample used to build a NICE extinction map. 
We used the T Tauri locus in a (H-Ks) v s. (J-H) colour-colour 
diagra m, defined by iMever et all (1 1 9971) and the ICardelli et al.l 
(Il989l) extinction law to remove all sources with (H-Ks) excess 
emission arising from the inner disk. 

The visual extinction Ay was then calculated for each star 
using the following equation 



15.93[(H -Ks) 



observed 



(H-Ks)i, 



cl 



(1) 



higher resolution (50") extinction map of the molecular cloud, 
shown in Figure |5] As for the 2MASS dust map, we removed 
disk sources from the photometric catalogue prior to building 
the FLAMINGOS dust map. As shown in section § 13.11 the re- 
gion covered by the FLAMINGOS data was also observed by 
Spitzer IRAC; we were therefore able to use the cirac disk iden- 
tification criterium (i.e., sources with ckirac > -2.56.) to remove 
these sources from the catalogue. 

To reduce contamination from foreground sources (that 
would erroneously lower the measured value of Ay for each 
pixel in the dust extinction map), we used only sources fainter 
than H=15 mag. The dust map of Figure |5] shows some "holes", 
in particular, one corresponding to the position of IRS 1 (~ 110 
southeast of the centre of the map); this hole is generated be- 
cause IRS 1 is a bright source that saturated at Ks and no good 
photometry could be extracted for the disk-less sources in its im- 
mediate vicinity. 

The densest regions in the FLAMINGOS extinction map 
correspond to the Spokes cluster (on which Figure |5] is cen- 
tred) and southwards towards IRS 1. The mean and median 
visual extinction values derived from the FLAMINGOS NIR 
dust map are 5.9 and 4.3 mag, respectively. Comparison of the 
FLAMINGOS NIR dust map with th e velocity integrate d '^CO 
(1-0) andC'**0 (1-0) maps obtained bv fRidge et alJ(l2003l) shows 
remarkable agreement with the structure of the molecular cloud. 
In particular, the molecular cloud component associated with the 
"Fox Fur Nebula" (a tilted "U"-shaped feature located in the 
northwestern part of the cloud) is clearly recognisable in the 
dust extinction map. The '^CO maps of Figure |6] shows that the 
cloud is essentially divided into two velocity components: (1) a 
4-6 km s"' component consisting of the bulk of cloud (including 
the Spokes cluster and the region surrounding IRS 1), and (2) a 
10-11 kms ' component corresponding to the Fox Fur Nebula. 

The mass of the molecular cloud was determined from 
the FLAMING OS dust extinction map, using the equation 
(lDickmanllI978l) 



where (H-K^) observed is the colour of a (disk-less) source ob- 
served towards the NGC 2264 field and (H-Ks)mirmsic is the in- 
trinsic stellar colour measured from a control field. Figure |4] 
shows the location of the control field used, which is represented 
by a box in the eastern part of the map. The area of the con- 
trol field is the same as that of the total area covered by the 
FLAMINGOS observations: 1265.48 arcmin^. To determine (H- 
Ks)intrinsic, wc uscd the median (H-K,) colour of stars in the con- 
trol field, i.e. 0.17 mag. The Ay of each star in a pixel (2.5' x 
2.5' box) was then averaged to determine the extinction for that 
pixel. The 2MASS dust extinction map has an average of 22 stars 
per pixel in the region of the molecular cloud. 

Since the 2MASS data is relatively shallow, it is unable to 
probe the densest parts of the cloud, which is why a pixel size 
of 2.5' X 2.5' was used: a finer resolution (smaller pixel size) 
yields a map with many "holes" corresponding to dense regions 
of the cloud where 2MASS did not detect background sources 
to the molecular cloud (i.e., empty pixels). As a result of the 
lower resolution of the dust map, the average Ay calculated for 
each pixel is an underestimate to the actual extinction because 
the higher Ay values are being diluted by those given by the 
neighboring less-extincted stars. 

Although the 2MASS extinction map is useful for an 
overview of the extent of the molecular cloud and surround- 
ing region, the resolution can be improved by a factor of three 
with the deeper NIR FLAMINGOS data. We therefore built a 



M = (ad)- 



,N(HI + H2) 



-^i 



2^v(0, 



(2) 



where a corresponds to the pixel size of the map in radians, d is 
the distance to the cloud, A^(///+//2)Mv=1.9 xlO"' cm"^ mag"', 
the mean molecular weight, fi, is 1.37, and the sum is 
over an individual pixel, /. The total mass obtained is 
M/oM/=3.7;^P gXlO^ Mq (the dominant contributor to the error in 
mass is the uncertainty in the distance). The mass contained in 
regions where Ay >7mag is M/i,,>7„,ai;=(2.1+O.5)xlO''M0, and 
the Fox Fur Nebula has Mfo^f^r^ (1.2+0.3)xlO^Mo. 

The total mass derived from the d ust extinction map^is con- 
sistent with the mass determined by iRidge et all (l2003i) from 
their '^CO (I-O) map, 5.2x10^ Mq (for a distance of 913pc), 
taking in account that their mass estimate is accurate to within a 
factor of 2- 4. Although the area covered by the FLAMINGOS 
dust map is 70 % that of the CO maps ( the '''CO map includes 
emission extending northwards and southwards of the dust map), 
our FLAMINGOS dust map covers all of the area with C'^O 
emission. Comparison of the mass in the dense regions of the 
cloud (MAy> 7maf;) with that dete rmined from the C'^O (1-0) map, 
2.1x1 O^Mm ( JRidge et al.l2003L for a distance of 9 1 3 pc). does in- 
deed show good agreement. 
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3.3. Extinction and spatial distribution of tiie disk population 

Figure Q shows two NIR+IRAC colour-colour diagrams of 
sources in the IRAC sample. The different source classes are 
clearly discriminated in these colour-colour diagrams: the disk- 
less sources (represented by small grey dots) occupy the region 
corresponding to the main sequence and its reddening band. The 
unembedded sources with disks are located along a band of 
relatively constant (J-H) colour Embedded sources with thick 
disks, flat spectrum sources, or Class I sources fan out in the 
NIR+IRAC diagrams at the red (KA^-S]) or (KA^]) end of 
the disk locus. These colour-colour diagrams may be used to 
identify and classify sources with disks (and circumstellar en- 
velopes). We note that the diagrams show a gap amongst the 
reddened sources, an effect that is particularly noticeable in 
panel (b) of Figure |7]- there are very few reddened sources with 
anaemic disks. 

To quantitatively determine the disk loci, we built four 
NIR+IRAC colour-colour density plots that are shown in Figure 
E (a) (/i:.r[3.6]), (b) (/:.r[4.5]), (c) (/:.r[5.8]), and (d) (/:,-[8]) 
vs. (J-H). We selected sources occupying a region in the (/^^-[S]) 
vs. {J-H) colour-colour space with a high density (>10mag"^) 
of excess emission sources to fit an empirical loci for T Tauri 
stars. The locus in a colour-colour diagram can be given by 



y — aQ + a\x, 



(3) 



where y corresponds to the {J-H) colour and x corresponds to 
one of the four colours: {Ks-i'i.b]), (ii:,-[4.5]), (/:,-[5.8]), or {K,- 
[8]). Table|3]summarises the coefficients aq and a\ obtained from 
the four fits. In the absence of spectral information, these loci 
may be used to identify embedded or unembedded sources with 
circumstellar disks. We compared the empirically determined lo- 
cus with the colour-colour distribution of unreddened models of 
disked sources in the Appendix A. Using a 5-cr limit for the disk 
locus of the {Ks-[^]) vs. {J-H) diagram, we are able to clearly 
separate sources with Ay greater than 3 mag. We note that pre- 
vious optical studies only identified PMS sources in the c luster 
for Ay <3 mag (e.g., see Figure 8 of lDahm & Simonll2005h . it is 
necessary to use longer wavele ngth data such as that presented in 
this p aper (and also shown in ISung et aDl2009l: ISung & Besselll 
I2OIOI) to identify the embedded population. 

Figure |9] shows the spatial distribution of sources with thick 
(a) and anaemic (b) disks relative to the dense material of the 
molecular cloud (represented by a Av= 8 mag grey contour). 
The thick and anaemic disk populations follow distinct spatial 
distributions, and also have distinct reddening distributions; the 
distribution of the embedded thick disk population is clearly 
segregated and c onfined to the densest regions of the cloud. 
ISung et aP (l2009l) 's analysis of the spatial distribution of sources 
with disks in NGC 2264 lead to the identification of three clus- 
terings: (i) near IRS 2 or the Spokes cluster, (ii) near IRS 1 or 
Allen's source, and (iii) near S Mon. For the sources with thick 
disks, groups (i) and (ii) are clearly identified in Figure |9] We 
additionally note that there is a loose clustering of embedded 
sources with thick disks in the Fox Fur Nebula region. The 
clustering is small compared to those of IRS 1 (Allen's source) 
and IRS 2 (Spoke's cluster) and is likely unbound; it is com- 
pose d of 19 sou r ces w ith thick disks, 9 of which are embed- 
ded. iRidge et aP (l2003h 's '^CO observations of the molecular 
cloud show that the bulk of the cloud has VLSR=4-6kms ' and 
that there is a separate component of Vlsr=10-11 kms"' cor- 
responding to the Fox Fur Nebula (see Figure |6). As this com- 
ponent collides with the bulk of the cloud, it could increase the 
local densities and subsequently trigger star formation, which 



could explain the clustering of embedded sources with the thick 
disks we identify here. 

Panel (b) of Figure |9] shows that most of these sources lie 
within the contour of the dense cloud material. Ay =8 mag, but 
otherwise there is no spatial segregation between sources with 
anaemic disks. Furthermore, sources with anaemic disks are not 
embedded, while sources with thick disks that are clustered (i.e., 
in the Spokes cluster, near Allen's source, and in the vicinity 
of the Fox Fur Nebula) are all embedded. ISung et al.l (120091) 
also defined a more dispersed halo population, more uniformly 
spread out in the cluster - its spatial distribution is consistent 
with that of the anaemic disk population. 

3.4. Accretion activity of the disk population 

To investigate dis k accretion, we consid er a subsample of Ha 
emitting sources dPahm & SimonI 120051) from our IRAC sam- 
ple; this subsample consists of 172 sources with measured Ha 
equivalent widths, W{Ha). Figure [TOl shows how W(Hq') varies 
for stars of different q-irac values; the figure is composed of 
four panels corresponding to stars of four spectral type groups. 
The W(Hq;) limit used to classify a sour ce as accreting or non - 
accreting is a function of spectral type (IWhite & Basril 120031) . 
These different limits are represented by vertical dashed lines in 
Figure [TOl sources located to the right of these lines are consid- 
ered to be accreting. 

The diagrams show that the majority (~85%) of the sources 
with thick disks appear to be accreting (consistent with what 
iLada et aD (l2006l) found in the young cluster IC 348), in compar- 
ison, 59% of sources with anaemic disks are accreting. Figure 
[TOl also shows the existence of sources with passive anaemic 
disks and sources with no IRAC disk that are non-accreting - 
these are expected to occur if the gas and dust removal in the 
disk occur at approximately the same rate. A particularly inter- 
esting region in the diagrams corresponds to that inhabited by 
accreting sources with no IRAC excess emission. There are 17 
sources with anaemic IRAC disks and 4 sources with no IRAC 
or inner disk that are accreting. Of the accreting sources with 
anaemic IRAC disks, six have excess 24yum detections indicat- 
ing the existence of thick outer disks, and fifteen of them have X- 
ray emission. Among the accreting sources with no IRAC disk, 
all four of them are X-ray emitting sources and three of them 
possess excess 24yum emission characteristic of a thick outer 
disk; the present data does not exclude the possibility that the 
other sources may also possess a thick outer disk. The X-ray 
sources with no IRAC disk and excess 24jum emission must 
therefore each have an inner hole that does not contain dust, 
but does contain some gas as it continues to accrete; the accre- 
tion is indicated by the large value of W{lia) and excess UV. 
We classify these types of sources as having transition disks, 
i.e., disks with inner dust holes filled with gas, and optically 
thick outer dust+gas disks. Accreting sources with anaemic in- 
ner disks and thick outer disks could be the immediate precursors 
to these transition disks. At this point, we remind the reader that 
the nomenclature for disk classification is still evolving and as 
such, non-standardized; a summary of the different terminolo- 
gies used by researchers can be found in [Evans et al. (2009). In 
particular, some authors would consider the disks that we de- 
fine in this paper as anaem ic to also be transition disks, because 
they are more evolved (e. g. Naiita et al1l2007l; ICurrie et alll2009l ; 



ICurrie & Sicilia-Aguilarll20Ill) . While some of the sources with 
anaemic IRAC disks that we identify in NGC 2264 may be pre- 
cursors to sources with disks that have cleared-out inner dust 
holes, others may not; we therefore decided to use a stricter def- 
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inition of the transition disk so as not to include anaemic IRAC 
disk sources. We return to this discussion in section § 14.21 

Well-known examples of sources with tr a nsition disks in- 
clude TW Hya (lOMvrl (ICalvet et alJ I2002L l20Q5h . GMAur 
(2-3Myr) (ICalvet et all 120051: iGrafe et al.l 1201 ll) . DH Tau. and 
DM Tau (e.g. lGrafe et al.ll201 lb . Transition disks have holes that 
are typically a few AU (< 10) in radius, and display a broad He? 
emission line. The spectroscopic study by Furesz et al . (2008) of 
sources in the Orion Nebula Cluster revealed sources that have 
no detectable disk emission short-wards of 8 //m yet show clear 
Ha accretion signatures. These authors found approximately 35 
sources with these characteristics, interpreting them as sources 
with transition disks. The age at which these sources may be 
found can therefore be as young as ~ 1 Myr 

Sources with no IRAC excess and Ha equivalent widths 
> lOA might also be dMe stars masquerading as sources with 
transition disks. Main-sequence M-type stars are classified as ei- 
ther having high (dMe) or low (dM) chromospheric and coronal 
activity; the magnetic activity is identified via Ha emission, i.e., 
dMe stars have typically equivalent widths 15>EW(Ha)>l A, 
while dM stars h ave EW(Ha)<lA (e.g. iHawlev et alj [19961 
lKruseetal1l2010h . In addition, dMe stars have Call HK reso- 
nance doublet lines in emission, as well as X-ray flaring. dMe 
stars may therefore be mis taken for e i ther p re-main sequence 
WTTS or transition disks. iRiaz et all (l2006l) analysed Spitzer 
IRAC and MIPS data for dMe stars and found that their SEDs 
could be described by a purely photospheric spectrum out to 
24//m, so longer wavelength data, i.e. >24//m, is necessary 
to ide ntify bona-fide transition disk sources. Dahm & SimojJ 
(l2005h discussed the contamination of dMe stars in their Ha sur- 
vey of NGC2264 and estimated it to be -8%; from their analy- 
sis, it is expected that foreground dMe stars are very faint and 
red in the (V-I) vs. I optical colour-magnitude diagram; after ex- 
amining optical photometry, we found that our sources do not 
populate this part of the diagram and are thus unlikely to be dMe 
stars. 



4. Discussion 

4.1. Star formation efficiency 

As previously mentioned, the FLAMINGOS NIR dust extinction 
map indicates that the molecular cloud associated with the clus- 
ter has mean and median Ay of 5.9 and 4.3 mag, respectively. 
We built a 2MASS control field iir^-band luminosity function 
(KLF) and reddened it by these two values, creating two red- 
dened KLFs, and subsequently subtracted each one of them from 
the FLAMINGOS KLF to statistically estimate the number of 
sources in excess that would thus be members of the NGC2264 
cluster We opted for calculating the population using both the 
mean and median values stated above, for these correspond to 
the lower and upper limits of the real value. 

Figure \n\ shows two KLFs ( the FLAMINGOS KLF and 
the 2MASS KLF uniformly reddened by Ay =4. 3 mag), as well 
as the corresponding differential KLF. The number of sources 
of the differential KLF is 1194, for A", <I5.5 mag (correspond- 
ing to a 0.2 Mq star, using lBaraff'e et al.l(ll998h PMS evolutional 
models, for a distance of 913 pc and an age of 2.5 Myr). Using 
the same procedure, but uniformly reddening the 2MASS KLF 
instead by Ay=5.9mag yields a population of 1677 sources 
(Ks <15.5mag) within the same area of 1264.48 arcmin^. The 
size of the stellar population of the NGC2264 cluster (within this 
area) is therefore estimated to be 1436+242. iLada et al.l ( Il993h 
used this same method to statistically determine the number 



of sources of NGC 2264; they found 360+130 sources {K,^ < 
13.5 mag), albeit for a smaller area: 540 arcmin^. 

Assuming, for s implicity, that each source has 0.5 Mq 
dMuench et al.l l2007h . then the stellar mass of NGC2264 is 
718+121 Mq. This value can be used to estimate the star 
formation efficiency (SEE) of NGC 2264, which is defined 
as SFE=stellar mass/(stellar mass + cloud mass). For the 
cloud mass, we used the mass measured from regions with 
extinction greater tha n the star formation extinction thresh - 
old, i.e.. Ay >7mag (iJohnstone et al.l 120041; iLada et alj|2010l) . 
(2.1±O.5)xlO^M0. The star formation efficiency of the cloud is 
calculated t o be ap proximately 25 %, similar to the value de- 
termined by iLada & Lada (2003) for nearby star forming clus- 
ters. A more accurate estimate of the total stellar mass can of 
course be obtained by using an appropriate IMF, although that 
is beyond the scope of this paper The reader should also note 
that the statistically determined members are all brighter than 
7*^5=15.5 mag, i.e. the fainter (deeply embedded young protostel- 
lar and/or lower mass) members are not accounted for and so the 
stellar mass used is an underestimate. 



4.2. Disk evolution 

A new empirical evolutionary paradigm is emerging for primor- 
dial optic ally thick disks, which consists of two distinct path- 
ways (e.g.lLada et al.ll2006l;lNaiita et al.ll2007HCieza et alJ | 2007 , 



l2008l JMerfti et all l2008t ICurrie et al.l l26o9riCieza et aiTl2010 
and references therein). In one of these evolutionary paths that 
a primordial disk can take, it experiences a h omologous deple- 
tion of dust (via dust settling and grain growth. lWood et al.l2002h 
throughout the entire extent of the disk; this path is characterized 
observationally by SEDs with slopes that do not vary substan- 
tially longwards of 3.6 yum. The alternative evolutionary path a 
primordial disk can undergo consists of a radial, inside-out clear- 
ing of the dust disk. The SEDs of sources that follow this latter 
path have different slopes shortwards and longwards of 8 fim. 
The results presented in this paper lend further support for this 
empirical disk evolution paradigm. 

We showed in section § l3.3l that the anaemic disk population 
differs from that of the thick disks in two main aspects: its spa- 
tial distribution displays no segregation, and it is an unembed- 
ded population (Ay <3 mag). Both of these facts indicate that 
the anaemic disk population is more evolved than the thick disk 
population, having had enough time to emerge from the parental 
molecular cloud and begin populating the halo of the NGC 2264 
cluster. 

Figure [12] characterizes the outer disk for the two afore- 
mentioned populations, using MIPS 24 /im data. The {Ks-[24]) 
colour is used as a proxy for estimatin g the optical thick ness 
of the outer disk (following Figure 2 of ICieza et al.l 120081) . We 
did not detect all of the IRAC sources at 24yum: only 29% of 
the anaemic sources have 24 //m counterparts, while 63% of the 
thick disk sources were detected at 24 //m. The 24 /im sources 
that are most likely detected are those with large excesses due 
to the cluster's distance, and indeed. Figure [12] shows that the 
IRAC disk sources detected at 24yum practically all have opti- 
cally thick outer disks. The diagram identifies sources with (i) 
optically thick inner and outer disks, (ii) anaemic inner disk and 
thick outer disk, (iii) anaemic inner disk and optically thin outer 
disk, (iv) sources with no inner disk and optically thick outer 
disks, and (v) sources with no inner disk and optically thin outer 
disks. 

The figure shows the path for a radial or inside-out evolu- 
tion of primordial optically thick disk, since we have identified 
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sources in all the different phases. The younger sources could 
evolve from (i) to (ii) and finally to transitional disks (iv). On 
the other hand, a homologous disk evolution could be described 
by young sources in (i), evolving into (iii). As explained before, 
region (iii) of the diagram is not well-populated by sources be- 
cause the 24fim detections are biased towards large excesses. 
Taking into account that 70% of the anaemic inner disks do not 
have lAfiva counterparts, we can reasonably assume that their 
outer disks are optically thin. These results may therefore indi- 
cate that the majority of primordial optically thick disks tend to 
follow homologous depletion evolutionary path. 

Multiple evolutionary paths have also been observationally 
identified in other young star-forming regions, suc h as IC 348 
(iLada et all l2006h . Lupus, Serpens, and Cha II (iMerin et al.l 
|2008|) . These regions are all much closer than NGC 2264, thus 
do no t suffer from the detection bias at 24 /^m. ( iLada et all 
l2006l) compiled an extensive atlas of SEDs for disk sources 
and found that the most numerous type of disk was one that 
can be described b y homologous dep letion. Furthermore, stud- 
ies of NGC 2362 fcurrie e t^l l2009h and the Coronet cluster 
dCurrie & Sicilia-A guilar 2 01ll) also show that the frequency of 
homologously depleted disks i s higher than that of disks with 
inner holes. On the other hand. lCieza et al.l (1201 Ol) carried out a 
detailed study of evolved disks in the nearby Ophiucus molecu- 
lar cloud, and argue that at one point during their evolution all 
disk sources pass through a phase where the outer disk remains 
optically thicker than the inner disk. 



4.2.1. Transition disks 

There are currently several possible scenarios that could ac- 
count for the formation of holes within protoplanetary disks. 
The first of these clearing processes is (aVphotoevaporation and 
visco u s disk evolution proces ses (e.g. [Alexander & Armitagg 
l2007h . [Alexander et al.l (l2006allbl) conducted a theoretical study 
of disk evolution, using a model combining both viscous 
evolution of the disk and photoevaporation of the disk by 
stellar radiation, and suggested that all evolving disks pass 
through a short inner hole phase; they furthermore predict 
that sources with transition disks should represent approxi- 
mately 10% of the disk sources. Another process that can 
clear an inner dust hole i s (b) grain growth and plan et for- 
mation (e.g. jj anaka et al.' '2005^. iHahertv & Muzerollel |2008|) . 
JDullemond & Dominik (2005) predict, using evolutionary mod- 
els, that grain growth and dust settling occurs more rapidly in the 
inner region of the disk than in the outer regions. Inner disk holes 
can also be formed through the (c) dyna mical interaction of the 
disk w ith an embedded giant planet (e.g. JAlexander & Armitaga 
|2009|) . Finally, the process of (d) dyna mical interaction of the 
disk with a (sub)stellar companion (e.g. Ilreland & Krausll2008h 
may also remove dust from the inner part of the disk. 

We identified four candidate transition disks from Figure 
[TOl these sources are accreting according to their Hff equivalent 
width and do not possess IRAC inner disks. We indicated three 
of these sources in Figure [12] with squares (the fourth candidate 
has no 24fim detection), and the diagram shows they all have 
optically thick outer disks. 

If we consider these four candidate transition disks, then the 
fraction of transition disks for NGC 2264 is ~4% (for a total 
105 disked sources with Ha measurements). Considering the 
24 //m detection bias in our data, our transition disk fraction 
may well be an underestimat e. This value is lowe r th an that 
found in other studies, such as ISicilia-Aguilar et al] (e.g. |2006|) 
who found that 10% of the disked sources in the young cluster 



Trumpler 37 are sources wi t h trans ition disks. This value is in ac- 
cord with lAlexanderet al.l (l2006al[bl) 's proposal that 10% of the 
dis ked sources are tran sition disks. A similar value was proposed 
bv lButleret 10(120061) . namely that -12% of FGKM stars host 
giant planets within 30 AU, implying that giant planet formation 
clears the inner disk. 



5. Summary 

We have found that sources with anaemic inner disks are more 
evolved than those with thick inner disks. We presented several 
results that lead to this conclusion, namely: 

1 . Sources with thick disks are mostly embedded, while sources 
with anaemic disks are mostly unembedded. Using the 
NIR+IRAC colour-colour diagrams and their disk loci, we 
determined the extinction towards sources with anaemic and 
thick disks: ~20% of the sources with thick disks are embed- 
ded (Ay > 3 mag), while none of the sources with anaemic 
disks are embedded (Ay > 3 mag). 

2. The spatial distribution of extincted sources with thick disks 
or sources with anaemic disks is very different: extincted 
sources with thick disks are found to group into three clus- 
terings (near IRS 1, the Spokes cluster, and in the Fox 
Fur Nebula), whereas sources with anaemic disks are more 
spread-out throughout the region with no clear clustering, 
populating the larger halo of the cluster This indicates that 
sources with anaemic disks have had time to move away 
from their birth-site and expand into the whole region. 

3. Most of the sources with thick disks have Ha emission that 
is consistent with accretion namely ~85%, while only ~59% 
of the sources with anaemic disks appear to be accreting. 

Our MIPS 24/Ym observations of the cluster detected 29% of 
the anaemic disk sources, and 63% of the thick disk sources. 
Unfortunately, our MIPS observations are very likely biased to- 
ward identifying sources with large excesses at 24yum due to the 
distance of the cluster. Most of the detected disk sources have 
optically thick outer disks. In addition, we identified four can- 
didate transition disk sources in NGC 2264. These sources are 
unembedded, appear to be accreting, and in three cases have op- 
tically thick outer disks. We also found that the sources with can- 
didate transition disks are more evolved than sources with thick 
disks. 

Our findings support the emerging disk evolution paradigm 
of two distinct evolutionary paths: one where the disk is homol- 
ogously depleted throughout its radial extent, and another where 
the disk evolves radially, from the inside-out. Since there are 
many more anaemic disks in NGC 2264 than transition disks, 
the homologous depletion disk evolutionary path appears to be 
the one that most primordial optically thick disks follow. 

Regarding the star formation history of the cluster, the spa- 
tial distribution of sources indicates that NGC 2264 may have 
undergone more than one star-forming event. The star forma- 
tion may have initiated in the north, near S Mon, and continued 
to progress toward the south where it is forming stars now in 
the Spokes cluster and in the vicinity of Allen's source. We also 
identified a potential triggering of star formation in the north- 
west of the cluster, in the Fox Fur Nebula. 

Finally, we calculated the mass of the dense molecular cloud 
associated with the cluster to be 2.1xlO''M0 (for Ay > 7 mag). 
We also estimated the total population of the cluster to be 
1436+242 stars, and a star formation efficiency of ~25%. 
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Fig. 1. Distribution of oirac for sources in the IRAC sample 
(open histog ram) and for 349 candidate members previously 
identified by iFuresz et al.l (l2006h (hatched-line histogram). The 
small inset is the same diagram but zoomed into the aiRAC re- 
gion corresponding to sources with thick disks. The shaded ar- 
eas correspond to the regions in the histograms occupied by 
Class I (C I) sources, flat spectrum (FS) sources, sources with 
thick disks (TD), sources with anaemic disks (AD), and sources 
with naked photospheres (NP). 
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Fig. 2. Comparison of V- and J-band luminosity functions for 
sources with thick disks (open histogram) and sources with 
anaemic disks (filled histogram). 
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Fig. 3. Variation in q-irac with spectral type. The horizontal thick 
lines correspond to the limits used to classify the evolutionary 
phases of protostars and PMS objects. 
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Fig. 4. Nyquist sampled, 2.5' resolution 2MASS dust extinction 
map of the molecular cloud associated with NGC 2264. The map 
is centred on {a, (5)(J2000.0)=(6Hl™02.6^ +09°36'12.2"). The 
dashed contour corresponds to Ay =4 mag, while the solid con- 
tour corresponds to Av=8mag. The box on the eastern part of 
the map indicates the size and location of the control field, while 
the larger box centred on the molecular cloud represents the area 
surveyed by IRAC. The FLAMINGOS fields, also centred on the 
cloud, are represented by boxes drawn with thicker lines. 
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Fig. 5. FLAMINGOS dust extinction map of the molecular cloud 
associated with the NGC 2264 cluster. The map is centred on (a, 
5)(J2000.0)=(6Hl'"02.6% +09°36'12.2"). The spatial resolution 
is 50"and the contours range from 6 to 30 mag in steps of 4 mag. 
For positional reference, the sources IRS 1, IRS 2, and S Mon are 
represented by filled circles, while the approximate southeastern 
edge of the Fox Fur Nebula is marked by a dashed curve. 
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Fig. 6. '''CO intensity maps, integrated for two velocity inter- 
vals, of the molecular clo ud associated with NGC 2264 clus- 
ter (see iRidge et aP l2003h . The bulk of the molecular cloud 
has Vlsr between 4 kms' and 6 kms', and a clear separate 
component is present at 1 0-11 kms' corresponding to the Fox 
Fur Nebula (compare with Fig. |5). The contours are the same 
for both panels, starting at L5Kkms"' with increments of 
L5Kkms-'. 
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Fig. 7. NIR+IRAC colour-colour diagrams of sources in the 
IRAC sample. The symbols are as follows: disk-less sources 
are represented by grey dots, sources with anaemic disks by 
black open circles, sources with thick disks by dark grey crosses, 
sources with flat spectrum SEDs by filled black stars, and fi- 
nally. Class I sources by open circles and a cross. The main se- 
quence is plotted as a black line. The diagrams show that sources 
with circumstellar material can be clearly distinguished from the 
main sequence population. Sources with anaemic disks populate 
the disk loci and are not reddened, generating a reddening gap 
that is particularly striking in panel (b). Sources with thick disks 
also populate the disk loci, but some are also reddened and fan 
out in the diagram along with flat-spectrum sources and Class I 
sources. 




K,-[X] 

Fig. 8. NIR+IRAC colour-colour Nyquist-sampled density di- 
agrams for sources in the IRAC sample. The beam size is 
0.3x0.3 mag and is represented by a pixel in the upper right cor- 
ner of panel (a). The contours are the same for each diagram and 
correspond to 3, 6.5, 10, 15, 20, 30, 50, 75, 100, 125, 150, 175, 
200, 250, 300, and 400 sources mag"^. The solid black curve 
corresponds to the main sequence of the Pleiades (see lLada et al.l 
2006). The dashed-dotted lines correspond to the loci of disked 
sources for each colour-colour space (see Table |3)- The dotted 
lines correspond to 3-cr (5-cr) error in the disk loci for the upper 
(lower) two panels. 
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Fig. 9. Spatial distribution of sources with thick disks (a) and 
anaemic disks (b). Embedded sources (A^ > 3 mag) are repre- 
sented by filled circles, and unembedded sources (Ay < 3 mag) 
are represented by open circles. The spatial extent of the molec- 
ular cloud is represented by a grey contour of Av=8 mag. Both 
maps are centred on IRS 2. For positional reference, the sources 
IRS 1, IRS 2, and S Mon are represented by plus signs, while the 
approximate southeastern edge of the Fox Fur Nebula is marked 
by a dashed curve in panel (b). 
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Fig. 10. Comparison of the Ha equivalent width (WHa) with 
ciRAc for sources of different spectral type (SpT). The horizon- 
tal lines mark the separation between disk-less sources, sources 
with anaemic, and sources with thick disks (see Table |2]i, while 
the vertical dashed lines separate accreting and non-accreting 
sources. 
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Fig. 11. Comparison of the NGC 2264 FLAMINGOS K,-hand 
luminosity function (KLF) (solid black line) to the reddened 
2MASS control field KLF (red hatched-line histogram). The 
histogram in the inset shows the differential KLF, with 1194 
sources. 



K,-[24] (meg) 

Fig. 12. (Ks-[24]) vs. ckirac diagram for sources with no IRAC 
disk (plus symbols), sources with IRAC anaemic inner disks 
(open circles), and sources with thick inner IRAC disks (filled 
circles). The plot is divided into nine sectors by grey dotted 
horiz ontal lines and vertical dashed lines (following ICieza et al.l 
l2008l) : see the discussion in § 14.21 for a detailed description of 
sectors (i) to (v). Sources that were identified as candidate tran- 
sition disks from Figure [TOj are additionally marked by open 
squares. 

Table 2. Source classification scheme using aiRAC- 



Source classification 



Q-iRAc value 



Class I (C I) sources 
Flat spectrum (PS) sources 
Sources with thick disks (TD) 
Sources with anaemic disks (AD) 
Sources with naked photospheres (NP) 



Q'IRAC> 0.5 

0.5 > Q',RAc> -0.5 

-0.5>Q'IRAC>-1.8 

-1.8 >Q'iRAc> -2.56 
Q'iRAc< -2.56 



Table 3. Coefficients for Equation |3] 



Colour-colour space 



ao 



ai 



(^,-[8]) vs. (J - H) 
(^,-[5.8]) vs. (y - //) 
(^,-[4.5]) vs. (J - H) 
(A:,-[3.6]) vs. (J - H) 



0.56+0.04 

0.47+0.03 

0.634±0.12 

0.40+0.13 



0.10+0.02 

0.25+0.03 

0.170+0.09 

0.65+0.10 
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Table 1. Summary of the FLAMINGOS observations. 
Field R.A. (J2000) Dec. (J2000) Passband Dithers Exp. time(s) Total exp. time (s) Observation date 



A-3 OeMl^ai.Sa^ +09°30'08.2" 



J 
H 



48 
14 
31 



60.0 
60.0 
35.0 



2304 27 Dec 2002, 31 Jan 2003, 01 Feb 2003 
840 31 Jan 2003 
1085 31 Jan 2003 



A-4 06''4r^31.33'^ +09°50'08.2" 



J 
H 



32 
32 
32 
30 



60.0 
60.0 
30.0 
35.0 



1920 27 Dec 2002, 01 Feb 2003 
27 Dec 2002, 01 Feb 2003 
2880 27 Dec 2002 
1050 01 Feb 2003 



~KA2 06'"40™ 10.22^ +09°30'08.2" 



J 
H 



14 
15 
14 
32 



60.0 
60.0 
60.0 
35.0 



840 
900 
840 
1120 



03 Jan 2003 
03 Jan 2003 

22 Mar 2003 

23 Mar 2003 



A-13 06''40" 10.22^ +09°50'08.2" 



J 
H 



16 
16 
16 

36 



60.0 
60.0 
60.0 
35.0 



960 
960 
960 

1260 



23 Mar 2003 
03 Jan 2003 

22 Mar 2003 

23 Mar 2003 



cl 



06''40"53.2P +09°51'08.0" 



J 
H 



32 
32 
66 

32 



60.0 
60.0 
30.0 
35.0 



1920 30 Jan 2003, 30 Jan 2004 
1920 28 Jan 2003, 30 Jan 2004 

28 Jan 2003 
3 100 29 Jan 2003, 30 Jan 2004 



c2 


06''4r'10.82'^ 


+09°31'06.8" 


J 


47 


60.0 


2820 


28 Jan 2003, 01 Jan 2005 








H 


16 


60.0 


960 


28 Jan 2003 








K. 


64 


30.0 


1920 


28 Jan 2003, 01 Jan 2005 
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Appendix A: Modelling the disk locus 



We discussed in section § l3.3l empiiical disk loci for NGC 2264. 
To better understand how the disk locus depends on stellar 
and disk par ameters, we used radia tion transfer models of disk 
sources from lRobitaille et aD(l2006l) to explore the (Ks — [8]) vs. 
(J - H) colour-colour diagram in more detail. Figure lA.ll uses 
models with disk inclinations <70° and varying stellar mass. 
The modelled sources are classified according to the SED slope 
into sources with either optically thick disks (black symbols) or 
anaemic disks (grey symbols), in the same manner as we did for 
our data (Table|2]i. For comparison, we overplot the empirically 
determined disk locus and 5-cr limits from Table |3]on each dia- 
gram. We found no significant variation in the locus for sources 
with disks with varying inclinations angles below 75°, although 
the locus does have significant dependence on the stellar mass in 
this colour-colour space. Figure lA!2l shows the same set of mod- 
elled sources as Fig. lA.ll but with inclination angles > 75° (i.e., 
closer to edge-on disks). These plots show us that the locus is 
not as well-constrained, and the sources fan out from the locus; 
this effect is stronger for sources with thick disks. 
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Fig. A.l. Color-color diagrams using models from 
[Robitaille et al. (2006) for disked sources with inclinations 
smaller than 75° for different ranges of stellar masses (Mq). 
The black solid line represents the Pleiades main sequence, 
the dot-dashed line corresponds to the disk locus, and the 
dotted lines represent the 5-cr error loci as determined by the 
observations of NGC 2264 (see §[33] Figure|8]and TableE). The 
black dots represent sources with optically thick disks, and the 
grey circles represent sources with anaemic disks. 
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Fig. A.2. Same as Fig. lA.ll except for models of disks with in- 
clination angles greater than 75°. 



the disk locus because it has a typical intrinsic width in {J-H) 
colour of ±0.3 mag (for a strong 5-cr cut), which translates to 
an Al/ of 2.7 mag. However, the diagram is useful for robustly 
identifying embedded sources and unembedded sources, as we 
demonstrated in section §13.31 




Fig.A.3. Synthetic color-color diagram built using the models 
show r i in Fig. lA.ll an d Fig.[A.2l with an IMF sampling dKroupal 
l200ll : lRobitaille et alJl2006l) . The solid line represents the main 
sequence, the dot-dashed line the empirical disk locus, and the 
dotted lines the 5-cr error loci discussed in § 13.31 and shown in 
FigurelDand Table 121 



To try and reproduce the empiral disk locus, we built 
synthetic cluster using the probabilities calculated by 



iRobitaille et alJ (l2006l i.e., assuming a Kroupa initial mass func 
tion between 0. 1 and 30 Mq, and with stellar ages sampled lin- 
early in time between lOOOyr and 2Myr, and random disk in- 
clination angles). Figure lAJJ shows the resulting {Ks - [8]) vs. 
(/ - H) colour-colour diagram. This plot shows that it is impos- 
sible to accurately determine extinctions of sources using solely 
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